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Hard magnetic ferrite with a gigantic coercivity 
and high frequency millimetre wave rotation 
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Magnetic ferrites such as Fe304 and Fe203 are extensively used in a range of applications because 
they are inexpensive and chemically stable. Here we show that rhodium-substituted £-Fe203, 
e-Rh;^Fe2-x03 nanomagnets prepared by a nanoscale chemical synthesis using mesoporous 
silica as a template, exhibit a huge coercive field (H^) of 27kOe at room temperature. 
Furthermore, a crystallographically oriented sample recorded an value of 31 kOe, which 
is the largest value among metal-oxide-based magnets and is comparable to those of rare- 
earth magnets. In addition, e-Rh;^Fe2-x03 shows high frequency millimetre wave absorption 
up to 209 GHz. e-Rho.uFe-] 86O3 exhibits a rotation of the polarization plane of the propagated 
millimetre wave at 220 GHz, which is one of the promising carrier frequencies (the window of 
air) for millimetre wave wireless communications. 
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Magnetic ferrite, which is composed of iron oxide, was 
discovered in the seventh century B.C. Since then, mag- 
netic ferrites have contributed to daily hfe^'^, for example, 
magnetic motors, magnetic recording media, magnetic fluids and 
electromagnetic wave filters. However, due to their low magneto - 
crystalline anisotropies^"^, the coercivity of magnetic ferrites is gen- 
erally low. Hence, the development of magnetic ferrites with a large 
coercive field (H^) is an important issue because this type of ferrite 
can be used as various advanced materials^"^. 

Of the several polymorphs of Fe203 (refs 9,10), Y-Fe203 and 
a-Fe203 phases have been well studied. Y-Fe203, which is a represen- 
tative soft magnet, is a ferrimagnet with a spinel structure. a-Fe203 
has a corundum structure and shows weak ferromagnetism due to 
the Dzyaloshinsky - Moriya mechanism^ In contrast, e-Fe203 is a 
rare phase, and research on this phase had been running behind^^" 
It is only recently that e-Fe203 was found to exhibit a huge coer- 
cive field of 20kOe at room temperature^^. Since then, its magnetic 
properties have received much attention ^^"^^ 

In this Article we synthesized rhodium- substituted e-iron 
oxide, e-Rh;^Fe2_;t03, nanomagnets and found that this series 
exhibits the highest value of 31kOe among metal oxides. 
This value is comparable to those of rare-earth magnets^^'^^. 
Furthermore, it displays electromagnetic wave absorption with 
the highest resonance frequency among whole magnetic materials. 
Herein, we report the synthesis, crystal structure, magnetic prop- 
erties, first-principles calculations, high frequency millimetre wave 
absorption properties and magnetic rotation of the propagated 
millimetre wave. 



Results 

Synthesis and crystal structure of 8-Rh^Fe2_;t03 nanomagnets. 

A series of rhodium- substituted epsilon-iron oxide nanoparticles 
(8-Rh^Fe2_;,03, x = 0, 0.04, 0.07, 0.11 and 0.14) were prepared 
by impregnating a methanol and water solution containing iron 
nitrate and rhodium nitrate into mesoporous silica, subsequently 
heating at 1,200 °C in air, and removing the silica matrix using an 
aqueous sodium hydroxide solution (Fig. la) (see Methods). The 
transmission electron microscope image indicates that the obtained 
sample is composed of nanoparticles with an average particle size of 
ca. 35 nm (inset of Fig. la; Supplementary Fig. SI). 

The X-ray diffraction patterns confirm that the crystal structure 
is orthorhombic with Pnali space group (Fig. lb). Rh substitutes for 
a specific iron site (Fe^ site) among the four nonequivalent iron sites 
(Fe^-Fei) sites) in the crystal structure of 8-Rh;^Fe2_;t03. The unit 
cell volume refined by the Rietveld analysis gradually increases with 
Rh concentration. At higher Rh concentrations, a-Rh;^Fe2_;c03 is 
formed as a byproduct (Supplementary Fig. S2 and Supplementary 
Table SI). The particle size distribution depends on the sintering 
temperature; the particle size increases as the sintering tempera- 
ture increases, and the phase changes from y-phase to e-phase and 
then to a-phase (Fig. Ic). The threshold sizes for the transformation 
from y- to 8-phase and from 8- to a-phase are evaluated to be 10 and 
55 nm, respectively (Supplementary Fig. S3). 

Magnetic properties. The magnetic properties were measured 
by a superconducting quantum interference device (SQUID) 
magnetometre. The temperature dependences of the magnetization 
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Figure 1 1 Synthesis, crystal structure and thermodynamic stability of 8-Rh;fFe2-x03 nanomagnets. (a) Upper row illustrates the synthesis of 
8-Rh^Fe2-x03. In the first step, a mesoporous silica template is synthesized in a conical flask to give a white powder. In the second step, an aqueous 
solution of Fe(N03)3 and Rh(N03)3 with methanol is impregnated into the silica to yield a yellow powder. In the third step, the resulting product is 
calcinated in air, and in the fourth step, the final product is obtained after etching the silica matrix with an aqueous NaOH solution. Lower left transmission 
electron microscope (TEM) image shows the obtained mesoporous silica where the small white dots are the mesopores. Illustration in the dotted frame 
depicts the formation process of e-phase during calcination; as the temperature increases, Fe(N03)3 and Rh(N03)3 in the mesopores turn into y-phase 
(-800 °C). Around the melting temperature of glass (1,000 °C), the particles begin to aggregate. As shown in the centre TEM image, e-phase is the main 
phase at 1,200 °C. Lower right TEM image is the final product after etching the silica matrix. The scale bars below the three TEM images indicate 30 nm. 
(b) (left) Crystal structure of y-phase. Green and light green indicate octahedral and tetrahedral Fe sites, respectively, (centre) Crystal structure of 
e-phase. Dark red, red, orange and yellow indicate the four nonequivalent Fe sites for Fe/\, Fes, Fe^ and Fe^, respectively, (right) Crystal structure of 
a-phase. Fe site is shown as blue octahedrons. Grey lines indicate the unit cell, (c) Calculated G//\/m,/ versus d curves for /-phases (/ = y, e and a). Green, 
red and blue curves correspond to y-phase, e-phase and a-phase, respectively. The curves are calculated under the condition of fly >iJi^> ii^, Gy <g^ < 
and {o^ - Oy)/{0(j^ - < (jd^- iJ^)/(id(^- ji^), where jij is the chemical potential of /-phase and cr,- is the surface energy of /-phase. 
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Figure 2 | Magnetic properties of 8-Rh;fFe2-x03 nanomagnets. (a) Magnetization versus external magnetic field plot of non-oriented e-Rho.MFe^ 35O3 
nanoparticles at 300 K exhibiting a coercive field of 27 kOe. (b) Magnetization versus external magnetic field plot of crystallographically oriented 
e-Rho.uFe^ 3503/resin at 300 K exhibiting a coercive field of 31 kOe. 



of E-^^¥e2-x^3 powder samples show that the Curie tempera- 
tures (Tc) are 505K (x = 0), 495K (x = 0.04), 484K {x = 0.07), 474K 
{x = 0.11) and 469 K {x = 0.14) (Supplementary Fig. S4a). The mag- 
netization versus external magnetic field plots at 300 K demonstrate 
that the value increases as the x value increases, that is, 22 kOe 
(x = 0), 23kOe (x = 0.04), 24kOe (x = 0.07), 25kOe {x = 0.n) and 
27kOe (x = 0.14) (Fig. 2a and Supplementary Fig. S4b). 

Furthermore, crystallographically oriented e-Rhg ^Fe 186^3 
nanoparticles dispersed in resin (see Methods) reveal a gigantic 

value of 31kOe (Fig. 2b), which is related to the magnetic ani- 
sotropy field (H^) as H^^HJ2 with a uniaxial magnetic anisotropy. 
The magnetic hysteresis loop indicates that the estimated value 
of the sample for x = 0.14 is 53 kOe, which is much larger than those 
of Y-Fe203 (0.2kOe) and Fe304 (O.SkOe). 

First-principles calculations. First-principles calculations of 
8-Rh;^Fe2_;t03 were conducted to understand the electronic struc- 
ture of this series (see Methods). The density of states (DOS) shows 
that e-Rh;^Fe2_;t03 is a charge-transfer type insulator with an 02p 
valence band and a Ve3d conduction band (Supplementary Fig. S5). 
Rh4d is distributed in the same energy region as the 02p band below 
the Fermi energy (£p) level, suggesting that Rh4d and 02p overlap 
greatly. In the occupied region, the 3d bands of the Feg and Fe^ sites 
are occupied by a spins, whereas those of the Fe^ and Fe^ sites are 
occupied by P spins, that is, (Fe^, Feg, Fe^, Fe^) = (j5, a, a, P) (Sup- 
plementary Fig. S6). Although we attempted several alternative spin 
configurations, their energies are higher than that of (P, a, a, P), sug- 
gesting that these other configurations are not suitable. This finding 
is consistent with the spin configuration reported for e-Fe203, which 
is a collinear ferrimagnet composed of positive sublattice magneti- 
zations at the Feg and Fe^ sites and negative sublattice magnetiza- 
tions at the Fe^ and Fe^ sites^^'^^. The charge density map shows 
strong hybridization between Fe^"^ and 0^~; the electron charge 
is partially transferred from to Fe^"^ through hybridization 
(Supplementary Fig. S7). 

Millimetre wave absorption properties. The electromagnetic wave 
absorption properties were measured under a zero-magnetic field at 
room temperature in the millimetre wave range (30-300 GHz) using 
terahertz time domain spectroscopy (THz-TDS) (see Methods and 
Supplementary Fig. S8). As shown in Fig. 3a, resonance frequencies 
(/r) are observed at 182 GHz (x = 0), 186GHz (x = 0.04), 190GHz 
(x = 0.07), 198 GHz (x = 0.n) and 209 GHz (x = 0. 14). These f, values 
are the highest among all magnetic materials. The/j. versus plot 



shows that the/j. value proportionally increases with increasing 
(the inset of Fig. 3a). Before this work, 8-Fe203 exhibited the high- 
est /j. value of 182 GHz^^. As for the quality of the absorption, the 
Q-factor of the absorption, which is expressed as Q=/j./A/, where 
A/ is full width at half maximum, is evaluated to be 38 (x = 0), 20 
(x = 0.04), 21 (x = 0.07), 19 {x = 0.n) and 17 {x = O.U). 

Magnetic rotation of the propagated millimetre wave. The rota- 
tion angle and ellipticity were measured under a zero -magnetic 
field using a magnetized 8-Rho.i4Fei 85O3 pellet sample, which was 
obtained by a strong pulsed magnetic field (ca. 8T). Irradiating with 
a millimetre wave from the N-pole side of the sample causes the rota- 
tion angle to display a dispersive spectrum centred at /^.^ 209 GHz 
where the minimum and maximum values are - 10° and +10° at 
196 and 220 GHz, respectively, whereas the ellipticity exhibits a peak 
at/j. with a value of + 0.3 (Fig. 3b, left). Upsetting the magnetic poles 
of the sample inverts the signs of the rotation angle and ellipticity 
(Fig. 3b, right). The observed phenomenon clearly originates from 
the remanent magnetization and its direction (Fig. 3c). The figure of 
merit (F), which is defined as F= (rotation angle (deg))/(absorption 
(dB)), has a high value of 20 degdB" ^ at 220 GHz. 

Discussion 

Usually, iron oxide (Fe203) transforms directly from y-phase into 
a-phase as the sintering temperature increases in the typical calci- 
nation process. However, in the present synthetic method, 8-phase 
exists as a stable phase due to the nanosize effects on the phase 
transformation via the following mechanism. In a nanoparticle, the 
free energy of each /-phase (G/, / = y, 8, or a) is expressed as the sum 
of the chemical potential (/i/) and surface energy (A/CT/), where Af 
and (7/ represent the molar surface area and the surface free energy 
of the particle, respectively^^; G/ = jU/ + A/CT/. Because Aj is equal to 
6V^j/d, where V^j and d represent the molar volume and particle 
diameter, respectively, the free energy per molar volume (G/ZV^i /) 
is expressed as G/ZV^i /^/i/ZV^i / + 6(7//ii. This equation suggests 
that the surface energy contribution increases as the particle dia- 
meter decreases. When the parameters in the equation satisfy the 
relationships oi jJi^ >jl^ > fl^, (Jy <a^ <(J^ and (cjg- (7y)/((7oc- G^) 
<{ji^-jjy)/{li(y^-li^), then the curves for Gy/Vj^y, Gg/Vj^g and 
Ga/^m,a intersect, leading to the stable 8-phase at an intermediate 
d value (Fig. Ic). Such a nanosize effect has also been reported in 
AI2O3 (refs 28,29) and Ti305 (ref. 30). 

Next, we considered the mechanism of the large value in the 
present material. One origin of the gigantic is that the particle 
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Figure 3 | Millimetre wave absorption and millimetre wave magnetic 
rotation of e-RhxFe2-x03 nanomagnets. (a) Millimetre wave absorption 
spectra of e-Rh^Fe2-x03 under a zero-magnetic field at room temperature 
for x = 0 (blue), x = 0.04 (light blue), x = 0.07 (green), x = 0.11 (orange) and 
x = 0.14 (red) measured by THz-TDS. Inset plots the resonance frequency 
(fr) versus coercive field (He). Colours of the dots are consistent with the x 
values of the absorption spectra. Black line is to guide the eye. (b) Millimetre 
wave magnetic rotation and ellipticity change of e-Rho.i4Fe-i geOs- Rotation 
angle (red dotted line) and ellipticity (black dotted line) spectra measured by 
irradiating from the N-pole side (left) and S-pole side (right) of the sample 
using THz-TDS. Solid lines are to guide the eye. (c) Schematic illustration 
of the magnetic rotation and ellipticity change of the propagated millimetre 
wave due to the optically induced magnetic dipole transition. The green and 
ochre arrows indicate the magnetic field and electric field components of the 
impulsive terahertz wave pulse, respectively. Red and blue sides of the disc 
denote the N- and S-poles of the magnetized pellet sample. 



of the 8-phase. Because shape anisotropy in the present series is 
negHgible (below 0.3 kOe), the observed value of 53kOe is 
determined by the magnetocrystalline anisotropy 

The results of the first-principles calculations indicate that the 
following mechanism explains the large magnetocrystalline anisot- 
ropy of this series. The charge density map (Supplementary Fig. S7) 
shows strong hybridization between Fe^ and ~ , that is, the elec- 
tron charge is partially transferred from ~ to Fe^ , resulting in 
an electron configuration ofd^^^ where q is the transferred charge. 
Such a hybridization causes a non-zero orbital angular momentum 
{L^O) on Fe ion, which induces a large magnetocrystalline anisot- 
ropy. Furthermore, Rh-O-Fe hybridization is observed around the 
Rh- substituted site, which enhances the magnetocrystalline anisot- 
ropy through the contribution from the orbital angular momentum 
on Rh. Thus, the large value of the present material is attributed 
to the nanoscale particle size, which is suitable to form a single mag- 
netic domain, and the huge value due to the magnetocrystalline 
anisotropy. 

Finally, we discuss millimetre wave absorption, magnetic rota- 
tion and ellipticity change of the propagated millimetre wave. In 
a ferromagnetic material with a uniaxial magnetic anisotropy, an 
electromagnetic wave is absorbed at/j., known as zero-field ferro- 
magnetic resonance (natural resonance). Because the magnitude of 
the value determines the/j- value, that is,fj. = vHJ2n, where v is 
the gyromagnetic constant, the observed high fj. values are due to 
the large values of e-Rh;^Fe2_;t03. In fact, the versus plot 
shows that the/j. value proportionally increases with increasing 
(Fig. 3a). Furthermore, zero-field ferromagnetic resonance gener- 
ates changes in the rotation angle and ellipticity of the propagated 
millimetre wave. When an impulsive terahertz wave pulse is irradi- 
ated into the sample, the magnetic field component of the terahertz 
wave tilts the spontaneous magnetization away from the magnetic 
easy- axis direction due to gyromagnetic effect, which causes the 
bulk magnetization to precess around the easy axis. This preces- 
sion produces a magnetic rotation of the propagated millimetre 
wave at a particular frequency (Fig. 3c; Supplementary Movie 1). 
These observed millimetre wave absorption, magnetic rotation and 
ellipticity change are due to the optically induced magnetic dipole 
transition, which differs from the conventional magneto-optic 
effect (Faraday effect) in the visible or near- infrared region through 
an optically induced electric dipole transition^^. 

In conclusion, we prepared a series of e-Rh;^Fe2_;t03 nanomag- 
nets via a nanoscale synthetic method using mesoporous silica as a 
template. e-Rho ^Fei 85O3 displays the largest value of 31kOe 
among metal-oxide magnets. The gigantic value has potential 
in future high-density magnetic recording media because the large 

value enables the particle size to be greatly reduced while main- 
taining coercivity. In addition, 8-Rh;^Fe2_;t03 exhibits a zero-field 
ferromagnetic resonance up to 209 GHz, and the magnetic rotation 
and ellipticity change of the propagated millimetre wave occur due 
to an optically induced magnetic dipole transition at 220 GHz. The 
present material should be useful for high frequency millimetre 
wave absorbers and rotators (isolators or circulators) as it should 
restrict electromagnetic interference problems^^ because the fre- 
quency of the magnetic rotation corresponds to the highest window 
of air (220 GHz band), which is the anticipated carrier frequency 
for next- generation millimetre wave wireless communications 
such as high -definition television communications and intra-board 
integrated circuit (IC) chip wireless communications^^"^^. 



size of ca. 35 nm in the present series should be sufficiently small 
to form a single magnetic domain^ '^^ In fact, the magnetic force 
microscopy image for x = 0 shows that each particle has a single 
magnetic domain (see Methods). In a single magnetic domain, only 
the magnetization rotation process operates, leading to a large 
value. Another origin of the gigantic is the magnetic anisotropy 



Methods 

Materials. An impregnation method using mesoporous silica nanoparticles 
was used to prepare a series of £-Rh;^Fe2_;c03 {x = 0, 0.04, 0.07, 0.11 and 0.14) 
nanoparticles. Mesoporous silica nanoparticles were synthesized by modifying 
reported methods^^'^^. Tetraethyl orthosilicate (62 mmol) was added to a solution 
of distilled water (4mol), ethanol (0.18 mol), cetyltrimethylammonium chloride 
(16 mmol) and triethanolamine (62 mmol) at 60 °C, yielding mesoporous silica 
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nanoparticles after centrifuging and heating at 600 °C in air. To obtain the 
target materials, a methanol and water solution containing iron nitrate and 
rhodium nitrate was immersed into mesoporous silica and sintered at 1,200 °C 
for 4 h in air. Silica matrices were etched by an aqueous sodium hydroxide 
solution at 70 °C. 

Physical property measurements. Elemental analyses of the prepared samples 
were conducted with inductively coupled plasma mass spectroscopy (Agilent 
Technologies, HP 4500). Transmission electron microscope images were 
acquired with JEM 2000EX. The 20-0 scan X-ray diffraction measurements 
were performed using Rigaku Ultima IV and Rigaku RINT2100 with Cu Ka 
radiation (A= 1.5418 A) at 293 K. Rietveld analyses were performed using the 
RIETAN-FP program^^. The magnetic properties were measured using a SQUID 
magnetometer (Quantum Design, MPMS 7). 

Preparation of an oriented sample. An oriented e-Rho.nFei 8603/resin 
was prepared. The E-Rhg ^Fei 85O3 nanoparticles were oriented by placing 
8-Rho i4Fei 85O3 nanoparticles dispersed ultraviolet curing resin (SD-2200, 
DIG Corporation) in an applied magnetic field of 12.7 kOe. The sample was then 
irradiated with ultraviolet light (EX250, Hoya Candeo Optronics Corporation) 
to harden the curing resin and to pin the oriented nanoparticles. 

Single magnetic domain of the nanoparticle. With respect to the magnetic 
domain, we previously performed magnetic force microscopy measurements on 
£-Fe203, which was prepared by a combination of reverse-micelle and sol-gel 
methods. The magnetic force microscopy image showed that £-Fe203 particles 
form a single magnetic domain. Therefore, £-'Rhx¥e2-x^3 considered to be 
composed of a single magnetic domain. Additionally, the magnetic domain sizes 
of common metal-oxide-based magnets are much larger than the size of our 
nanoparticles. For example, the magnetic domain size of barium ferrite is 
about 1 [im^^. 

First-principles calculations. First-principles calculations were performed using 
Vienna ab initio Simulation Package, which is a plane-wave projector augmented 
wave method program'^ The exchange-correlation functional was approximated 
with generalized gradient approximation (GO A) parameterized by Perdew-Burke- 
Ernzerhof (PBE) in the spin-restricted condition. To describe the ¥e3d and Rh4d 
orbitals, the Hubbard (7 term (GGA+ (7) was considered. Referring to the reported 
calculations for a-Fe203 and 'Y-Fe203, the difference between the Coulomb 
interaction energy (U^-f^^i) and exchange interaction energy (/gx), ^;oul~/ex' 
set to 4.0 eV. The basis set size was regulated by a cutoff energy of 520 eV. Gaussian 
methods (smearing width: 0.2 eV) were used to integrate the Brillouin zone with 
a 3x3x3 /c-mesh and to calculate the DOS with a 5x3x3 /c-mesh. Among the 
16 Fe atoms in the unit cell, a Rh atom substituted for one Fe atom at the Fe^ site. 
Supplementary Figures S5-S7 show the calculated DOS and charge density map. 

Millimetre wave absorption measurements. A pump-probe measurement 
system was used for the THz-TDS measurements. Supplementary Figure S8 
schematically diagrams the THz-TDS system. A mode-locked Ti:sapphire pulse 
laser delivered 20-fs long ultrashort light pulses with a central photon energy of 
1.55 eV at a pulse repetition rate of 76 MHz. A beam splitter divided the output of 
the mode-locked Ti:sapphire femtosecond pulse laser into pump and probe beams. 
Dipole type and bowtie type low- temperature-grown GaAs photoconductive 
antennas were used as the THz wave emitter and detector, respectively. The emitted 
THz wave was coUimated and irradiated into the sample, which was mounted on 
a sample holder and placed between a set of paraboloidal mirrors. The transmitted 
THz electric field was detected with the bowtie type low-temperature-grown 
GaAs photoconductive antenna, which was triggered by optical probe pulses with 
variable time delays between the pump and probe pulses; consequently, the electric 
field of the transmitted THz pulse wave was formed in the time domain. The 
absorption spectra of the THz waves were calculated as (Absorption) = - 10 log 
|t(fo)p [dB], where t{(o) is the complex amplitude transmittance. Samples were 
held in paper containers with filling factors of ca. 30%. 

Millimetre wave magnetic rotation measurements. To obtain the millimetre 
wave magnetic rotation spectra, we employed three polarizers (P1-P3) to measure 
the horizontal and vertical polarization components (Supplementary Fig. S8). The 
measurement system operated only the horizontal polarization component of the 
electric field because the antenna determined the directions of the THz wave emit- 
ter and detector. To precisely define the polarization of the THz wave as horizontal 
(0°), PI was placed after the emitter while P3 was placed before the detector. The 
components of + 45° and - 45° were measured by rotating P2 behind the sample. 
The vertical and horizontal polarization components were obtained as their differ- 
ence and sum, respectively. The rotation angle and ellipticity were obtained from 
the horizontal and vertical components {E^{t), Ey{t)) of the transmitted electric 
fields, which are expressed as 

(t) = jai (co) cos{cot + d{co))dco, Ey (t) = j ^2 (^) cos(ft)f + 82 {co))dco. 



The rotation angle (p{(o) and ellipticity r}{(o) spectra are expressed as 



2 yaf{co)-a2{co) 



77(^0) = tan isin-M ^^^^^^^sind{o,) \\, 



where 5{co) = di{(o) - d2{co) is the phase difference. a{co) and 5{(o) can be directly 
determined as the Fourier-transformed amplitude and phase of the temporal 
waveforms. 
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